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Abstract 

Many binary minor planets (BMPs; both binary asteroids and binary Trans-Neptunians objects; 
TNOs) are known to exist in the Solar system. The currently observed orbital and physical properties 
of BMPs hold essential information and clues about their origin, their evolution and the conditions 
under-which they evolved. Here we study the orbital properties of BMPs with currently known 
mutual orbits We find that BMPs are typically highly inclined relative to their orbit around the sun, 
with a distribution consistent with an isotropic distribution. BMPs not affected by tidal forces are 
found to have high eccentricities with non-thermal eccentricity distribution peaking at intermediate 
eccentricities (typically 0.4 — 0.6). The high inclinations and eccentricities of the BMPs suggest that 
BMPs evolved in a dense collisional environment, in which gravitational encounters in addition to 
tidal and secular Kozai affects played an important role in their orbital evolution. 



1. INTRODUCTION 

The binary asteroids and binary Trans-Neptunian ob- 
jects (TNOs) discovered in recent years show a large 
diversity of orbital properties, showing a wide range of 
eccentricities, inclinations, separations and mass ratios. 
Many models have been suggested for the origin of these 
binary minor planets (BMPs) and their orbital configu- 
rations (Rich ardson fc W alsh 2006). An essential com- 
ponent in constraining theoretical models for the origin 
and evolution of BMPs is understanding the distribu- 
tion of their orbital parameters. The number of BMPs 
with known orbital parameters is currently small. Nev- 
ertheless, 29 BMP systems already have full solutions for 
their mutual orbits, including 17 TNOs and 12 asteroids 
(we do not consider near earth objects which have much 
shorter lifetimes). We study the orbital properties of the 
29 main belt and transneptunian binaries (which we re- 
fer to corporately as BMPs) to provide clues to and con- 
straints on their evolutionary history. Several reviews 
have presented the o bserv ed separati ons of BMP s (e.g . 
IRichardson fc Walshl[2006l : INoll et, al.l[20M IWalsrJl2009h 
here we focus on the distributions of eccentricities and 
inclinations, not shown before. In addition, we discuss 
the relations between the orbital parameters of BMPs, 
including their observed periods/separations. 

In the following we present and discuss the published 
orbital solutions of BMPs. We present the distributions 
of the orbital parameters of these BMPs, treating binary 
TNOs and binary asteroids separately. We also discuss 
the selection biases affecting both samples. We then 
briefly study the implications of our findings regarding 
the conditions in the early Solar system and the forma- 
tion and evolutionary scenarios of BMPs. 

2. THE ORBITAL PARAMETERS OF BMPS 
2.1. The data 
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The mutual orbits of 30 BMPs have been published in 
the literature (see tables 1 and 2), some of them with 
two degenerate solutions, and a few published with no 
indicated inclination. Table 1 shows the physical prop- 
erties and orbital parameters of BMPs with known incli- 
nations in the solar system. In the literature the orbital 
parameters of BMPs are typically given with respect to 
the ecliptic plane, while the inner orbital parameters are 
given with respect to the equatorial frame of reference. 
For our analysis we are interested in the mutual inclina- 
tions between the BMP orbit around the Sun and the 
BMP inner orbit, therefore we used the published data 
to calculate these inclinations for our analysis of the dy- 
namical evolution (see Appendix). 

In some cases two degenerate orbital solutions were 
found for the BMPs; these solutions differ significantly 
in their derived inclinations, but have very similar ec- 
centricities. In these cases we detail both solutions (see 
table 1) and we use the published eccentricities to de- 
rive the eccentricity distribution (shown in Fig. [I}. All 
published binary asteroids orbits have unique solution 
for the inclination (besides the binary asteroids Balam 
for which only the eccentricity is known) and we show 
their mutual inclination distribution (in cosi; Fig. [2]). 
Only 5 of the binary TNOs, however, have a unique non- 
degenerate solutions for their inclinations. We therefore 
can not show the true inclination distribution for binary 
TNOs with significant statistics (Fig. [2] shows the distri- 
bution of binary TNOs inclination for an arbitrary choice 
of one solution for each binary TNO from the two pos- 
sible degenerate solutions published). Nevertheless, we 
do consider all the possible distributions of the binary 
TNOs in a statistical manner (see section 3). Note that 
both binary TNO and binary asteroid populations are 
presented. Given that these two populations differed in 
the conditions under which they evolved, their orbital 
properties are presented separately. 

2.2. Selection effects 

There are significant observational selection effects 
present in the current distribution of binary TNOs and 
asteroids. These effects are very difficult to correct since 
they depend on several factors, and the currently ob- 
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Fig. 1. — Eccentricity distribution of observed binary asteroids 
and TNOs (panels (a) and (b) respectively). Panel (c) shows the 
cumulative distribution of both samples (asteroids and TNO bina- 
ries, solid lines, respectively) as well as comparison to a thermal 
distribution (lower solid line). The dashed line shows the cumu- 
lative eccentricity distribution of binary TNOs excluding binaries 
that are likely affected by tides (e < 0.05). 
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Fig. 2. — Inclination distribution of observed binary asteroids 
and binary TNOs (panels (a) and (b) respectively). Panel (c) shows 
the cumulative distribution of both samples as well as comparison 
to a flat distribution (i.e. random in cosi). Note that many of the 
binary TNO inclinations have two degenerate solutions; here only 
one solution is chosen arbitrarily for each binary to illustrate their 
likely distribution (see text). 

served population of binaries comes from very heteroge- 
neous observing conditions. Although we do not correct 
for these effects in the current analysis, we discuss them 
below to recognize the possible biases they may produce. 

2.2.1. Binary TNOs 

Binary TNOs are primarily discovered in two ways. 
Ground-based observations can detect two co-moving 
objects, if the projected separation of the objects is 
roughl y larger than the t ypical seeing, >0.5 arcseconds 
(e.g. iKern fc E lliot 2006) o r if adaptive optics is used on 
bright targets (e.g. iBrown et al.ll2006D. The wides t bi- 
naries, such as 2001 QW322 (|Petit fc Mousisl 120041) can 
be found with this mechanism, assuming that discov- 
ery and follow-up images have been searched for bina- 
ries (generally true, but not always). Pan-STARRS and 
other future s urveys will probably detect dozens of bina- 
ries this way (|Holman et al.lr2007l i. The second discovery 
mechanism is using Hubble Space Telescope (HST) ob- 
servations combined with PSF-fitting techniques, which 
has discovered the great majority (~ 80 %) of known bi- 



nary TNOs. HST can resolve binaries as close as tens 
of milliarcseconds and generally pro bes much fainter ob- 
jects than ground-based telescopes (|Noll et al.ll2008l) . 

Each of these two methods suffer from important ob- 
servational biases. The main bias is a detection bias: the 
secondary must be bright enough and far enough away 
from the primary to be detected. In both the ground- and 
space-based cases, there is no simple prescription that 
describes this detection bias, especially since these are 
coupled close to the primary: a bright secondary can be 
discovered at smaller separations than a faint secondary. 
Furthermore, since HST observations are always nearly 
the same duration (one HST orbit), observations are es- 
sentially magnitude-limited, implying that low bright- 
ness ratios (and mass ratios) can only be seen around 
brighter targets. 

A discussion of these biases and t heir effects on th e 
observed population is also given bv INoll et all ((20081 ). 
Despite the observational biases, these authors believe 
that there is good evidence that binaries composed of 
moderate-sized TNOs (diameters less than ^1000 km) 
are actually clustered at nearly equal brightness, with 
A„ lag < 1 corresponding to a mass ratio of greater than 
^0.25 assuming equal albedos and densities. This lat- 
ter assumption seems reasonable given that the colors 
of components of bina ry TNOs are known to be similar 
(jBenecchi et al.ll2009aft . 

Considering now biases in the binary mutual orbital el- 
ements, we again point out that the majority of binaries 
are discovered in single HST snapshots. Hence, biases are 
introduced by the fact that the orbital separ ation must 
be det ectably large at a single epoch. While N oll et al.l 
(2008) point out that the statistical distribution of ob- 
served separations is similar to the distribution of semi- 
major axes, it is clear that systems with larger eccentric- 
ities are more likely to be seen since the observability is 
increased at apoapse, both due to wider separations and 
to the longer residence time. This eccentricity bias is 
most important near the angular resolution limit where 
most of the binaries are discovered, but it is not im- 
portant when semi-major axes exceed ~0.2 arcseconds. 
Of course, separation-limited observations imply that the 
smallest semi-major axes are undetectable; these systems 
are best discovered through photometry, either through 
doubly periodic light curves or where eclipses and oc- 
cultations may reveal con tact binaries like 2001 QG298 
(jSheppard fc JewittlfeOO^ . 

Finally, we also note a bias against observations of bi- 
nary TNOs with large size/mass ratios, due to detection 
limits. For example, a typical TNO with a radius of 100 
km located at 40 AU would have collisional satellites with 
V < 24.5 at separations of ^0.04 arcseconds, perhaps 
barely detectable by HST if the satellite is at elonga- 
tion. The most likely way to discover these systems is 
through photocenter-barycenter shifts (in this case, the 
size would be ~2 milliarcseconds) detected through long- 
baseline highly-accurate astrometry, potentially available 
from future Pan-STARRS or LSST surveys. Detecting 
a double-periodic light curve (or serendipitous mutual 
events) may be possible for some systems, though the 
faintncss of the components makes this very difficult and 
the result may be impossible to distinguish from a single 
object with arbitrary shape and spin orientation. The 
size ratio of the binary components can be important for 



Orbital properties of binary minor planets 



3 



the origin and evolution of the binary orbital parame- 
ters, and therefore may affect the statistics of the orbital 
parameters in our sample. 

Another important bias is that systems with low in- 
clination^ present nearly edge-on orbits with respect to 
Earth-based observations, while systems with high rela- 
tive inclinations are usually seen face-on. Since binaries 
are discovered when the components are significantly sep- 
arated on-the-sky, there is a greater likelihood for low 
inclination secondaries to be unresolvable when observed 
at a single random epoch. This bias is reduced as the 
projected semi-major axis grows, but remains significant 
even at a few times the resolution limit. 

2.2.2. Binary asteroids 

While some binary asteroids have been imaged by 
HST, ground-based adaptive optics has been employed 
more often for these brighter systems than their trans- 
Neptunian counterparts. The much smaller sizes of typ- 
ical asteroids is partly offset by their increased bright- 
ness and proximity. Discovery of binaries photometri- 
cally through double-periodic light curves a nd /or mutual 
event s is common for near-Earth asteroids (jPravec et al.1 
120061 note, however that these are not included in 
our analysis, and are mentioned here for completeness). 
Other methods are radar observations that often reveal 
near-Earth binary asteroids (not discussed here), and 
stellar occultations that ca n reveal the presence o f main 
belt binary asteroids (e.g. iDescamps et al.ll2007[ ). The- 
oretically, the radar and occultations techniques have 
fewer observational biases than the other more common 
techniques, but their application is severely limited. 

Binary asteroids are also clearly subject to detection 
bias: objects with smaller satellites are more difficult to 
observe. However, since asteroids are searched for bina- 
ries using ground-based facilities, there is a greater pos- 
sibility of searching for companions in more than a sin- 
gle snapshot. Furthermore, the geocentric orientations 
of these systems change much more rapidly than for es- 
sentially fixed KBO orbits (~ 60° yr -1 for asteroids, vs. 
~ 2° yr^ 1 for KBOs). Therefore, the eccentricity and 
inclination biases are not as strong as in the Kuiper belt. 
In binaries discovered through mutual events, there is an 
obvious bias towards edge-on systems, though "edge-on" 
can probe a wide range of inclinations. 

3. DISCUSSION 

Several different processes affect the mutual orbits of 
BMPs. Some are related to their initial formation and 
others to their later evolution either as isolated systems, 
or due to the effects of external perturbations and en- 
counters with other objects. The studies of these pro- 
cesses have been focused on a specific type of BMPs such 
as binary asteroids or binary TNOs. However, all of these 
suggested processes could in principle be relevant for the 
formation/evolution of BMPs both close (binary aster- 
oids) and far (binary TNOs) from the sun. 

The various suggested mecha nisms for the f orma- 
tion of BMPs (see refs. lAstakhov et all 120051 : 

4 Here, inclination means the mutual inclination between the 
heliocentric orbit and the mutual binary orbit; for objects at these 
great distances, the difference between heliocentric and geocentric 
viewing angles is not significant for this bias. 



Richardson & Walshl [20061 : INoll et all 120081 for some 
overviews) predict different initial orbital configurations. 
These include the following: 

• Smashed t a rget satellites (S MATS ; 

iWeidenschilimel [20021: iDurda et al.l [2004 [20ll : 
Low eccentricity distribution expected from 
collisionally formed satellites orbiting the main 
collision remnant body (which typically form as 
close binaries and are likely to be affected by 
tides). 

• Escap e ejecta binaries (EEBs; IDurda et al.l 120041 . 
2010): Intermediate eccentricities from bound 
ejecta pairs ejected following a collision of two 
larger bodies. 

• Exchanged binaries (jFunato et al. 2004) : Typically 
high eccentricities (> 0.8) for high mass ratio bi- 
naries formed through exchanges. 

• Dynamical fricti on and chaos assisted capture 
(CAC ) binaries (|Goldreich et all 1200% iLee et all 
120071 ): Intermediate eccentricities (0.2 < e < 0.8) 
for binaries formed through chaos assisted capture 
of satellites. 

Unfortunately, only a few studies expl ored aspects of 
the inclination distribution o f BMPs (lAstakhov et aTl 
2001 iNazzario et all [20071 iSchlichting fe Saril 120081 
Perets k, Naozll2009j ). 

After their formation, BMPs can be affected by several 
processes which can change their orbits. Tidal effect are 
most important when the BMPs components approach 
each other at a close distance. These effects couple the 
orbital evolution of the BMP to the spin of the BMPs 
components, and the tides raised on the objects serve 
to dissipate the total angular momentum of the system. 
Tidal effects can also excite and enla rge the eccentricity 
and i nclination of a given BMP (e.g. iGoldreich fc Soterl 
1966). At long enough timescales (depending on spe- 
cific configuration), however, tidally evolved systems are 
expected to relax into more circularized configuration, 
possibly locked configurations and even mergers. Such 
effects are thought to produce the period-eccentricity dis- 
tributions of (close) binary systems such as stellar bi- 
naries and planetary system s, and are lik ely to play a 
similar role in BMP systems ([Mazeh 20Q8|). 

Another evolutionary process is the Kozai-Lidov 
ijKozail Il962t iLidovl I1962D mechanism, which is the ef- 
fect of a secular perturbation from a third object (in a 
triple system, i.e. the Sun serves as the third compan- 
ion for BMPs) on the (bound) binary system. It could 
lead to large (order unity) periodic oscillations (Kozai 
cycles) in the eccentricity and inclination, i.e. it could 
both raise and lower the inclinations and eccentricities 
of a system. Note, however, that such a process is ef- 
fective only for systems with initially high inclinations 
((40° > i m < 140°; with a somewhat wider inclination 
range for initially eccentric systems). 

The combined effects of the Kozai mechanism in ad- 
dition to tidal friction, (Kozai cycles and tidal fri ction: 
KCTF lMaieh fc ShaharrJ[l979l:lKiseleva et alJll998l) . can 
change the orbital parameters of the BMPs, and re- 
duce both the eccentricity and separation of the BMPs 



TABLE 1 
The orbital parameters of BMPs. 



Satellite Period 
Name (days) 



Separation 
(kn 



i(J2000) 
(deg) 



D rati* 



(deg) (deg) 



(kn 



Q 
(deg) 



Ref. 



3749 Balam* 




61 ±10 


289 ± 13 


0.90 








0.43 






15,25 


45 Eugenia 


Pet it- Prince 


4.77 ±0.001 


1180 ± 8 





109 ± 2 


126.78 


112 


0.81 


202 


203 ± 2 


16,17 


22 Kalliope 


Linus 


3.6 ±0.001 


1095 ± 11 





99.6 ± 0.5 


103.42 


-92.5 ± 60 


0.15 


181 ±4.6 


284.5 ± 2 


16,17 


283 Emma 




3.35 ±0.00093 


581 ± 3.6 


0.12 ± 0.01 


94.2 ± 0.4 


65.46 


40 ± 4 


0.06 


160 


345.4 ± 0.4 


16,25 


130 Elektra 




5.26 ±0.0053 


1318 ± 24 


0.13 ± 0.03 


25 ± 2 


23.75 


311 ± 5 


0.04 


215 


1.6 ± 2 


15,25 


379 Huenna 




87.6 ±0.026 


3335.8 ± 54.9 


0.22 ± 0.01 


152.7 ± 0.3 


168.84 


284 ± 5 


0.06 


92 


204.3 ± 0.3 


15,25 


762 Pulcova 




4.44 ±0.001 


703 ± 13 


0.03 ± 0.01 


132 ± 2 


131.9 


-189 ± 20 


0.14 


137 ±3.2 


235 ± 2 


16,25 


90 Antiope 




0.69 ±4.1 x 10~ 6 


171 ± 13 


0.03 


63.70 ± 2 


54.74 


60 ± 30 


0.95 


87.8 ±1 


303.1 ± 2 


6,17 


121 Hcrmione 




2.56 ±0.0021 


747 ± 13 





79.1 ± 4 


70.19 


84.3 


0.17 


187 ±68 


83.7 ± 3 


7,13 


107 Camilla 




3.72 ±0.003 


1250 ± 10 





17 ± 5 


28.32 


-32 


0.06 


249 


141.00 ± 2 


15,25 


87 Sylvia* 


Romulus 


3.65 ±0.0007 


1356 ± 5 





7 


27.89 


-87 ± 11 


0.46 


282 ±4 


101 


13,25 


617 Patroclus 


Menoetius 


4.28 ±0.004 


680 ± 5 


0.02 ± 0.02 








0.92 


60.9 




17 


42355 Typhon 


Echidna 


18.97 ±0.0064 


1628 ± 5 


0.53 ± 0.02 


37.9 ± 2 


50.56 


09 


0.55 


76 t\t 


253.1 ± 4 


10,21 


19990J4 




84.09 ±0.016 


3303 ± 5 


0.37 ± 0.01 


53.80 ± 1.2 


119.56 


53.96 




37.5 ±8.5 


275.8 ± 2.2 


11 






84.14 ±0.016 


3225 ± 18 


0.36 ± 0.01 


99.8 ± 1.5 


56.7 


71.7 




37.5 ±8.5 


210.2 ± 1.6 


11 


90482 Orcus 


Vanth 


9.54 ±0.0001 


8980 ± 18 





90.2 ± 0.6 


92.13 





0.31 


900 


50 ± 0.6 


4 






9.5392 ±0.0001 


8985 ± 24 





305.8 ± 0.6 


70.02 





0.31 


900 


249.4 ± 0.4 


4 


47171 TC36* 




1.9068 ±0.0001 


867 ± 11 


0.101 ± 0.006 


88.9 ± 0.6 


74.21 


77.7 


0.93 


265±« 


330 ± 1 


1 


134340 Pluto* 


Charon 


6.39 ±10 -6 


19571.4 ± 24 





96.16 


119.61 





0.49 


2302 


223.05 ±10~ 4 


3,19,21,25 


134860 2000OJ67 




22.04 ±0.004 


2361 ± 36 


0.09 ± 0.02 


84.6 ± 3 


85.21 


-233.9 




69 ±16 


272.9 ± 3.1 


11 






22.04 ±0.0036 


2352 ± 35 


0.09 ± 0.02 


73.80 ± 2.9 


94.45 


136.8 




69 ±16 


212.2 ± 3.3 


11 


2001XR254 




125.61 ±0.12 


9326 ± 75 


0.56 


41.07 ± 0.22 


20.27 


-94.76 




84.5 ±19.5 


341.16 ± 0.33 


11 






125.61 ±0.13 


9211 ± 69 


0.55 


154.50 ± 0.22 


155.36 


-21.88 




84.5 ±19.5 


125.18 ± 0.55 


11 


136108 Haumca* 


Hi'iaka 


49.13 ± 0.03 


49500 ± 69 


0.05 


234.8 ± 0.4 


139.14 


278.6 ± 0.4 






26.1 ± 0.4 


2,21,24 


66652 Borasisi 




46.26 +S'™5 


4660 ± 170 


0.46 ± 0.01 


152 ± 3 




159.94 ± 4.01 


1 


316 




20,25 






46.23 

— 0.074 


4700 ± 170 


0.45 ± 0.01 


51 ± 3 




167.39 ± 3.44 


1 


316 




20,25 


2001QW322 




9855 
6570 


114000 ± 
105000 ± 


0.2 
0.4 


118 
130 








54 
58 




22 
23 


88611 Toharonhiawako 


Sawiskcra 


876 ±227 


31409 ± 2500 


0.31 ± 0.08 


128.1 ±6.5 


128.8 


330.3 ± 22.4 


0.69 


78 


96.70 ± 13.4 


22,25 


2003TJ58 




137.32 ±0.19 


3799 ± 54 


0.53 ± 0.01 


38.1 ± 2.1 


62.25 


-110.04 




32.5 ±7.5 


194.60 ± 4.2 


11 






137.32 ±0.19 


3728 ± 44 


0.53 ± 0.01 


96.1 ± 2 


116.77 


-88.90 




32.5 ±7.5 


150.80 ± 2.8 


11 


1998WW31 




574 ±10 


22300 ± 44 


0.82 ± 0.05 


41.7 ± 0.7 


51.96 


159.50 


0.83 


118 


94.30 ± 0.8 


25,26 


2004PB108 




97.02 ±0.07 


10400 ± 130 


0.44 ± 0.01 


89 ± 1.1 


84.13 


229.93 




120.5 ±27.5 


121.99 ± 0.75 


11 






97.08 ±0.069 


10550 ± 130 


0.45 ± 0.01 


106.55 ± 0.99 


95.23 


211.91 




120.5 ±27.5 


30.19 ± 0.86 


11 


58534 Logos 


Zoe 


312 ±3 


8010 ± 80 


0.45 ± 0.03 


121.5 ± 2 




310.13 ± 2.87 


0.825 


80 




19,25 






310 ±3 


7970 ± 80 


0.37 ± 0.01 


69 ± 2 




298.09 ± 5.73 


0.825 


80 




19,25 


2000QL251 




56.46 ±0.018 


4991 ± 17 


0.49 ± 0.01 


127.78 ± 0.62 


135.7 


42.2 




74 ±17 


109.5 ± 1.1 


11 






56.44 ±0.017 


5014 ± 16 


0.49 ± 0.01 


45.62 ± 0.66 


46.58 


45.70 




74 ±17 


71.20 ± 1.1 


11 


136199 Eris 


Dysnomia 


15.77 ±0.002 


37430 ± 140 


0.01 


61.3 ± 0.7 


94.98 






2400 ±100 


139 ± 1 


3,8,21 






15.77 ±0.002 


37370 ± 150 


0.01 


142 ± 3 


85.74 






2400 ±101 


68 ± 3 


3,8,21 


65489 Ceto 


Phorcys 


9.55 ±0.007 


1841 ± 47 


0.02 


116.6 ± 3 


115.38 


-64.6 


0.17 


87 


134.6 ± 3.4 


9,23 






9.56 ±0.008 


1840 ± 47 


0.01 


68.8 ± 2.9 


66.28 


-65.5 


0.79 


87 


105.5 ± 3.7 


9,23 



o 

bfj 
-. 

=3 



CD 



o 



References: 1 IBenecchi etTai! 120091)1) 2 IBrown et aTJ 120051) 3 IBrown fc S challcr ( 201171) 4 IBrown et ail 12010D 5 IBuie et all 12006D 6 IDescamps~ 
8 IGreenberg fc Barnei 120081) 9 fGrundv et"al] 12007D 10 IGrundv et al.l 120081) 11 IGrundv et al.l 120091) 12 IHestroffer et al.l 120051) 13 IMarchis et al.l 
14 IMarchis et al.l 120061) 15 IMarchis et al.l J2008H) 15 IMarchis et al.l 12008a!) 17 IMargot fc Brownl 120011) 18 INoll et aTJ 120031) 19 INoll et al] 12004a! 
21 INoll et al.l 120081) 22 IQsip et all 120031) 2 1Petit et all 12008TI 24 [Rabinowit z et al.l 1200g) 25 IRichardson &: Wal sh ( 2006) 26 IVeillet et aTJ (2002) 



aLl 120071) 7 |pcscamps et al. (2009) 
(2005) 

INoll et al.l 12004bfl 



Errors estimates are also shown where available. Unnamed satellites are omitted. Multiple systems are marked with (*) and the orbital parameters are for the (listed) outer satellite. 
All measured orbital parameters are given with respect to J2000 equatorial plane. Calculation method of the mutual inclinations is given in the appendix. 
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TABLE 2 

The external (helicentric) orbital parameters of BMPs 



Name 


SMA 
(AU) 




1-out 

(deg) 


(deg) 


Mass 
(10 18 kg) 


error 
(10 18 kg) 


class 


rcf. 


3749 Balam* 


2.24 


0. 


.11 


5.39 


295.84 





±0.00002 


FF 




45 Eugenia* 


2.72 


0. 


.08 


6.61 


147.92 


5.69 


±0.12 


MB 




22 Kalliope 


2.91 


0. 


.1 


13.71 


66.23 


8.10 


±0.2 


MB 




283 Emma 


3.04 





.15 


8 


304.42 


1.38 


±0.03 


EF 




130 Elektra 


3.12 





.21 


22.87 


145.46 


6.6 


±0.4 


MB 




379 Huerma 


3.13 


0. 


.19 


1.67 


172.07 


0.38 


±0.019 


TP 




762 Pulcova 


3.15 


0. 


.1 


13.09 


305.8 


1.4 


±0.1 


MB 




90 Antiope 


3.16 


0. 


.16 


2.22 


70.22 


0.83 


±0.02 


TF 




121 Hermione 


3.44 


0. 


.14 


7.6 


73.18 


4.7 


±0.2 


OMB 




107 Camilla 


3.48 


0. 


.08 


10.05 


173.12 


11.2 


±0.3 


OMB 




87 Sylvia* 


3.49 


0. 


.08 


10.86 


73.31 


14.87 


±0.06 


OMB 




617 Patroclus 


5.22 





.11 


22.05 


44.35 


1.36 


±0.11 


JT 




42355 Typhon 


37 


.65 





.53 


2.43 


351.96 


0.95 


±0.052 


Cent 




19990J4 


38. 


.10 


0. 


.02 


2.61 


127.46 


0.40 


±0.0087 


ICC 


11 


90482 Orcus 


39. 


.17 


0. 


.23 


20.58 


268.65 


632 


±5 


3:2N 


4 


47171 TC36* 


39. 


.7 


0. 


.23 


8.41 


97.08 


14.2 


±0.05 


3:2N 




134340 Pluto* 


39 


.45 





.25 


17.09 


110.38 


14570 


±9 


3:2N 




134860 2000OJ67 


42 


.9 





.01 


1.33 


96.76 


2.15 


±0.099 


cc 


11 


2001 XR254 


43 







.02 


2.66 


52.73 


3.92 


±0.089 


cc 


11 


136108 Haumea* 


43. 


.08 





.2 


28.22 


122.1 


4200 


±100 


HF 




66652 Borasisi 


44 


.07 


0. 


.09 


0.56 


84.74 


3.8 


±0.4 


CC 




2001QW322 


44. 


.28 





.02 


4.8 


124.67 


1.5 


+0.3 
-0.2 


CC 




88611 Teharonhiawako 


44. 


.29 





.02 


2.57 


304.63 


3.2 


CC 




2003TJ58 


44 


.5 





.09 


1.31 


37.12 


0.22 


±0.0078 


CC 


11 


1998WW31 


44. 


.64 





.09 


6.81 


237.1 


2.7 




CC 




2004PB108 


45. 


.1 


0. 


.11 


19.19 


147.38 


9.88 


±0.37 


HC 


10 


58534 Logos 


45. 


.5 


0. 


.12 


2.9 


132.64 


0.42 


±0.02 


CC 




2000QL251 


47. 


.8 





.21 


5.83 


223.29 


3.14 


±0.03 


2:1N 


10 


136199 Eris 


67 


.96 


0. 


.41 


43.97 


35.99 


16600 


±200 


SD 




65489 Ceto 


100.17 





.82 


22.32 


172.04 


5.42 


±0.42 


Cent 





Unless noted otherwise all outer parameters are taken from JPL small bodies database (see http://ssd.jpl.nasa.gov/sbdb.cgi). Errors are 
shown where available. Multiple systems are marked with *. References are detailed in Table 1. Ine heliocentric orbital classification 
is noted as "class", where we used the following notations: MB, OMB, JT, CC, ICC, HC, SD, Cen, HF, TF, EF, FF and n : mN for 
Main Belt, Outer Main Belt, Jupiter Trojan, Cold Classical, Inner Cold Classical, Hot Classical, Scattered Disk, Centaur, Haumea Family, 
Themis Family, Eos Family, Flora Family for n : m Neptune resonances respectably. Asteroid family membership is based on listings in 
I Zappala et a l. (1995). The prevalence of binaries among the cold classical population of the transneptunian belt is discussed further in 
INoII et alj (2008)" 
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(jPerets fe Naodl2009l ). In essence this effect rapidly low- 
ers the eccentricity of BMPs and shortens their period. 

These mechanisms are effective in isolated systems 
(although including the Sun). In collisional systems, 
encounters between BMPs and other minor planets or 
BMPs can change the orbital parameters of the BMPs. 
The distribution of inclinations in such systems is likely 
to be randomized, where as the eccentricity distribution 
is expected to appr oach high eccentricities, on average 
(|Funato et alJ 12004ft, poss ibly producing a thermal like 
distribution (Heggiei ri975ft . 

In the following we discuss the implications of the ob- 
served distributions in the light of the dynamical pro- 
cesses discussed above. 

3.1. Eccentricities 

The observed eccentricity distribution of BMPs (Figs 
Q]) shows both low (and zero) eccentricity BMPs as well 
as high eccentricity ones (up to 0.82 and 0.22 for binary 
TNOs and asteroids, respectively). The clear correla- 
tion between the eccentricity and the semi-major axis of 
the BMPs (eccentricity-period distribution of BMPs is 
shown in fig [3]), reminiscent of o ther binary p opulations 
(e.g. binary stars, exoplanets; iMazehl [20081 ) indicates 
that low eccentricity BMPs are likely to be produced by 
tidal circularization, which become important for BMPs 
of small separations. In principle, the period eccentric- 
ity distribution can be used to constrain tidal evolution 
theories and/or the physical parameters of BMPs which 
affect the tidal evolution (e.g. the Q parameter, and its 
evolution Efroimsk v fc Lainevl [2007) . Although current 
statistics are still too small to produce strong constraints 
on such theories/parameters, one can already check spe- 
cific tidal evolution cases (see e. g. the theoretical lines 
shown in fig. [3J. 

The apparent lack of zero and low (< 0.2) eccentric- 
ity binary TNOs at larger separations (where tidal forces 
are not effective) suggests that either the formation pro- 
cesses of binary TNOs are not inclined to form them at 
such eccentricities, e.g. in the EEB and CAC scenarios 
for BMPs formation, and/or that later dynamical evo- 
lution changed their eccentricities. The current sample 
of binary TNOs, showing the lack of many high eccen- 
tricity BMPs (e > 0.8), is already large enough to rule 
out the exchang e formation s c enario for binary TNOs 
as formulated bv lFunato et al.l ((2004) as the main single 
process producing their current distribution. We note 
that the lack of high eccentricity for the largest binary 
TNOs together with the frequent intermediate eccentric- 
ity population of the lower mass binary TNOs could 
be consistent with a collisional scenario. In this pro- 
cess EEBs are made from smaller remnant bodies follow- 
ing a collision where as the bigger remnants may form 
SMA Ts with lower eccen tricities in non-disruptive coll- 
sions (|Durda et al.ll2004f) . 

The eccentricity distribution of binary asteroids ap- 
pears to differ from that of binary TNOs, with a larger 
fraction of binaries at shorter more circular orbits. How- 
ever the general trend is similar to that of binary TNOs, 
showing close binaries to typically have more circular or- 
bits and wider ones to have higher eccentricities. Such 
distribution could be consistent with that of binary 
TNOs, given the selection effects (e.g. the difficulty in 
finding close small sized binary TNOs) and the small 



statistics. Moreover, the lack of wider period binaries 
with higher eccentricities is likely related to the much 
smaller phase space available to binary asteroids, due to 
their smalle Hill radii (see fig. 3). The binary aster- 
oids distribution may therefore also be suggestive of a 
collisional origin. The known collisional families in the 
asteroid belt give further support for this scenario. 

3.2. Inclinations 

The distribution of BMPs inclinations shows a large 
fraction of them to have high inclinations. The inclina- 
tions of binary asteroids are consistent with a random 
distribution of inclinations (flat in cosi). The underly- 
ing true inclination distribution of the larger sample of 
binary TNOs can not be derived directly (given the de- 
generate inclinations solutions for most of the sample). 
Nevertheless, we can statistically verify whether it too 
could be consistent with a flat distribution. To do so, 
we consider all the possible inclination distributions of 
the binary TNOa (i.e. 2 Nd '<>, where N deg = 9 is the 
number of orbits with two degenerate solutions used). 
We then use the Two-sample Kolmogorov-Smirnov test 
to check whether each possible distribution is consistent 
with it being drawn from a flat distribution (in cos i). We 
find that all of these distributions are consistent with an 
isotropic distribution. We conclude that the inclinations 
of both binary TNOs and asteroids are consistent with a 
random distribution of inclinations (flat in cosi); clearly 
BMPs ar e not restricted to plana r co nfigurations as sug- 
gested bv lGoldreich et al.l (12002ft and ISchlichting fe Saril 
( 2008). In some scenarios (|Goldreich et al.1120021) . BMPs 
are formed in a thin planetesimals disk with low velocity 
dispersion. It was suggested that BMPs with high in- 
clinations are therefore not likely to be produced under 
these conditions. However, three body encounters can 
easily change the inclinations of BMPs, as these could be 
highly chaotic, and produce highly inclined orbits even 
under such conditions (Perets & Kupi, in prep.). Col- 
lisoinaly formed BMPs are also likely to produce a range 
of inclination, as material could be ejected in a wide 
fan, and specifcially EEB could have random inclinations 
(however further studies in this direction are required). 
The high inclinations of BMPs could therefore be sugges- 
tive of the collisional environment at which BMPs were 
formed. Note that since encounters between BMPs and 
other single planetesimal can easily erase the initial dis- 
tribution of BMPs inclinations, predic tions of the inclina- 
tion d istribution such as suggested by ISchlichting fc Saril 
( 2008) , which do not seem to be consistent with the cur- 
rently observed distribution of inclinations (i.e. non pla- 
nar configurations) , are not likely to constrain formation 
scenarios of BMPs. 

We note that the combined processes of secular Kozai 
evolution (due to perturbations by the sur0) and tidal 
friction (KCTF) can lead to specific correlations between 
inclinations and separations of BMPs (as well as eccen- 
tricities). These processes and the currently observed re- 
lations they may have produced, hav e been discussed in 
detail elsewhere (|Perets fc N aoz 2009, see also Ragozzine 
& Brown, in prep.). We find that the probability for the 

5 Note that alternatively/in addition planets or even a third 
companion in triple minor planets systems could also produce such 
perturbations. 
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eccentricity and inclinations (|cosi|) distributions to be 
uncorrelated is ~ 0.04 (with the correlation coefficient 
found to be 0.7). 

We conclude that although understanding the shape 
of the inclinations distribution requires more data, it 
is clear that high inclinations serve as the rule and 
not the exception, and point to an excitation mecha- 
nism, possibly a highly collisional environment at the 
epoch of BMPs formation and/or evolution, with fur- 
ther KCTF evolution playing a role and producing cor- 
relations between the BMPs orbital properties (separa- 
tions /eccentricities / inclinations) . 

3.3. Dependence on mass/size 

The masses and sizes of BMP components, mostly the 
smaller secondaries, are not accurately known. Never- 
theless, some interesting trends with masses/sizes may 
already be observed in our binary TNOs sample. Binary 
TNOs with the largest primaries (> 500 km; see table 1 
and fig. 3) , are observed to have satellites at small sepa- 
rations and low eccentricities, possibly indicating a colli- 
sional formation mechanis m (|Canuplf2005t iBrownl [20081 : 
iRagozzine fc B rown 200jl). We also note that these bi- 
naries seem to cluster at relatively high inclinations (not 
shown, but see table 1), although this may not be sta- 
tistically significant given our currently small sample, it 
may suggest KCTF (highly efficient at high inclinations) 
was involved in catalyzing mergers /collisions or tidal evo - 
lution of pre-formed binaries (see iPerets &: Naoz| [2009). 
We do not see any trends in the current sample of binary 
asteroids, but note that no massive (> 500 km) binary 
asteroids are known. Nevertheless, it is worth pointing 
out that most of the binary asteroids share similarities 
with the massive binary TNOs, i.e. small separations 
and low eccentricities, as discussed above, suggesting a 
similar, possibly collisional, origin. 

3.4. Binary trans-Neptunian objects vs. binary asteroids 

We find that both the observed populations of binary 
asteroids and binary TNOs present similar orbital prop- 
erties (e.g. high inclinations). This suggest that some ba- 
sic features in their formation and evolution were similar 
(e.g. dense collisional environment). Although binary 
TNOs show typically much higher eccentricities, most 
of these binaries have much wider orbits than those of 
binary asteroids (both due to observational selection ef- 
fects, as well as the much smaller Hill radii in which bi- 
nary asteroids can exist; see figure 3). The differences in 
eccentricities may therefore only reflect the tendency of 
closer binaries to be more circularized (since tidal friction 
effects become stronger ). 

3.5. Multiple systems 

The last few years have seen the discovery and char- 
acterization of the first asteroids and TNOs with mul- 
tiple satellites, which deserve special mention. Such 
systems with well-kn own published orbits include Pluto 
( Thol en et al.1 l2008h. Haumea jlRagozzine fc Brownl 
I2009D. 1999 T C36 IBenecchi et al.l 12009b!) and Svlvia 
(|Marchis et al.l l2006fT while the asteroids Eugenia, 
Balam, Kleopatra, Minerva, 2001 SN263, and 1994 CC 
have only recently been announced as triple systems. 
Though the detection biases described in Section 2 are 




log (a/r p ) 



Fig. 3. — Eccentricity as a function of the binary separations 
normalized to the primary's radius. We consider both TNOs (open 
circles) and asteroids (squares). The large circles correspond to the 
largest size primary TNOs (> 500 fern). The solid line shows the 
critical separation-eccentricity beyond which typical binaries are 
expected to be strongly affected by tides and be circularized (where 
we choose Pulcova and OJ4 as representative for the asteroid and 
TNO binaries respectively) . The dotted lines show the the critical 
separation-eccentricity beyond which typical binaries become un- 
stable due to the perturbation by the Sun (the Hill radius); lines 
are shown for both asteroids (left line) and TNOs (right), as the 
Hill radius is strongly dependent on the distance from the Sun.. 

present for these systems, the majority of these systems 
were first known as binaries with additional companions 
found during subsequent study. It is therefore difficult to 
estimate the frequency of multiple systems in the various 
populations. Nevertheless, the existence of several such 
objects indicates that these are not rare. Detailed obser- 
vations of these systems can yield mass determinations 
for each of the bodies independently, which is not possi- 
ble for binaries (|Tholen et all 120081: IRagozzine fc Brownl 
120091) . 

Multiple systems provide unique additional leverage 
in determining the formation and evolution of binaries. 
For example, the coplanar nature of the satellite sys- 
tem of Pluto likely requ ires a dense collisional formation 
([Stephens fc Nolll 12006). although the detailed forma- 
tion and evolution of this system is still not understood 
(IWard & Canud 120061 : lUthwick fc Wull200l . Even for 
systems with unknown orbits, the small sizes and com- 
pact configurations suggest that all of these systems are 
collisionally formed, except for 1999 TC36; though mul- 
tiple episodes of YORP -induced fission m ay be relevant 
for the smaller bodies ( W alsh et al.l 120081 ). The hierar- 
chical and nearly-equal mass nature of the 1999 TC36 
triple system cannot be explained by a single collision 
and is likely the result of sequential formatio n by capture 
(e.g. iGoldreich et al.ll2002l : fLee et al.ll2007h . Balam has 
a very unusual satellite system, including an outer satel- 
lite wi th a putative eccentricity of ^0.9 (Marchis et al. 
2008b) , and possibly an unbound satellite t hat separated 
from Balam less than a million years ago (jVokrouhlickvl 
I2009T ). 

These systems also present examples of unique orbita l 
evolution. As pointed out by (Ragozzine & Brown 2009), 
the combination of rapid orbital expansion (compared 
to the weak expansion around giant planets, which are 
ineffective at dissipating tidal energy) and gravitation- 
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ally interacting satellites, creates a unique brand of tidal 
evolution. Multiple resonance crossings can excite ec- 
centricities aiid_2ncliir^iojis J _perhaps leading to instabil- 
ity (e.g., iCanup et al.l 11999 1. In hierarchical systems 
with significant mutual inclinatio n, the Kozai-Lidov ef - 
fect may destabilize the system (jPerets fc Naoz| [2009). 
For all multiple systems, survival to the present epoch 
can be difficult, and multiples must have been more com- 
mon in the primordial population. 

Additional study of these multiple minor planet sys- 
tems will provide unique insight s into the forma t ion and 
evolution of these systems. See iBenecchi et al.l {2009b) 
for an additional review of multiple minor planet sys- 
tems, including a table of properties. 

4. SUMMARY 

In this paper we compiled a catalog of BMPS, both 
binary TNOs and binary asteroids, with full orbital so- 
lutions. We presented a first analysis of the eccentricity 
and inclinations distribution of BMPs as well as their 
semi-major axis-eccentricity distribution. This data and 
its analysis can be used to study and constrain forma- 
tion and evolutionary scenario of BMPs. Specifically we 
find high relative inclinations for the BMPs as well as 
typically large (but not extremely high) eccentricities for 
BMPs not affected by tidal evolution. By themselves 
these results already suggest BMPs evolved in dense en- 
vironment in which collisions and close gravitational en- 



counters formed/perturbed the binaries and strongly af- 
fected their orbital evolution. We suggest that these en- 
counters together with secular Kozai evolution and tidal 
effects, could have erased much of the direct signatures 
of the initial formation of BMPs, as possibly reflecte d 
by their observed orbits (see also IPerets fc Naozl 12009). 
More theoretical work, however, is required to under- 
stand the different parts played by the initial formation 
and configurations of BMPs vs. their later dynamical 
evolution. Especially important are better theoretical 
predictions for the observational signatures of different 
BMP formation scenarios, which could be compared with 
the data presented in this paper and additional future 
data. 
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APPENDIX 

CALCULATION OF THE MUTUAL INCLINATION 

The mutual inclination i m represent the angle between the inner and the outer orbit. All of the heliocentric 
parameters are given with respect to the ecliptic plane, while the inner orbital parameters are given with respect to 
the equatorial frame of reference. Thus, we first transform the heliocentric inclination and the longitude of ascending 
node from the ecliptic plane (i c ,f2 c ) to the equatorial frame of reference (i q ,£l q ). Let us define the pole vector with 
respect to the ecliptic (equatorial) frame, P c (P 9 ): 













-( 









sinri C]g sinj Ci g 
— cos £l c ,q sini Cj( j 



(Al) 



The rotation matrix from the ecliptic to the equatorial is rotation with respect to the x axis of the invariable plane, 
i.e., 



1 







R(x) c ^q = ( cose —sine 
sin e cos e 



(A2) 



where e is the ecliptic angle. Thus, the resulting equations of the transformation from the ecliptic to the equatorial 
are: 



cos(i q ) = cos(i c ) cos(e) — sin(z c ) sin(e) cos(fi c ) , 
cos(Q q ) sin(i g ) = cos(i c ) sin(e) + sin(i c ) cos(e) cos(r2 c ) , 
sin(O g ) sin(iq) = sin(i7 c ) sin(i c ) , 



(A3) 



where from the last two equations we can find Q q using the atan2(x,y) function. 

Now we turn to calculate the mutual inclination. We take the longitude of ascending node and inclinations with 
respect to equatorial frame (resulting from the above transformation) of the inner and outer orbit, i.e., i out , riout, and 
i s ,£l s , and calculate the mutual inclinations. Again let us define the binary outer orbit and inner pole vectors as: 



sin f2 ut,in sin i ut,in 

COS ^out,in Sin 2 ut,in 
COS 2out,in 



(A4) 
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and thus, the mutual inclination is simply given by 

i m = cos -1 (P out • P in ) . (A5) 
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